Ballistic and diffusive: theory of vortices in the two-band superconductor MgB 2 
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Motivated by the recent results on impurity effects in MgE>2 , we present a theoretical model for a 
two-band superconductor in which the character of quasiparticle motion is ballistic in one band and 
diffusive in the other. We apply our model to calculate the electronic structure in the vicinity of 
an isolated vortex. We assume that superconductivity in the diffusive (n) band is induced by that 
in the clean (a) band, as suggested by experimental evidence for MgB2. We focus our attention to 
the spatial variations of the order parameter, the current density, and the vortex core spectrum in 
the two bands. Our results indicate that the coupling to the 7r band can lead to the appearance of 
additional bound states near the gap edge in the a band that are absent in the single-band case. 

PACS numbers: 74.20.-z, 74.50.+r, 74.70.Ad, 74.81.-g 
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It is now well established that MgB2 is a two-gap su- 
perconductor, and its essential superconducting proper- 
ties are well described byan isotropic s-wave two-band 
model Q, S H S S H B- The two 'bands' in MgB 2 
are the 'strong' a band that arises from the boron a- 
orbitals (with the energy gap « 7.2meV), and the 
'weak' 7r band that derives from the boron 7r orbitals 
(Ajr » 2.3meV). Despite the fact that each 'band' con- 
sists in reality of a pair of bands, the variation of the gap 
within each corresponding pair of Fermi-surface sheets 
can be neglected 0, . The two energy gaps are observed 
to vanish at a common transition temperature T c : no sec- 
ond transition has been observed, and there is evidence 
of induced superconductivity in the it band |a, llOl . 

There have been considerable efforts to understand im- 
purity effects in MgB2 0] ■ Besides the potential applica- 
tions of MgB2 as magnetic devices 01 > these studies have 
aimed at testing a prediction of the two-band model: the 
reduction of T c and the gap ratio by non-magnetic im- 
purities. However, Mazin et al. jl'lll 1| have shown that 
this does not apply to MgB2, i.e. to Mg-site impurities 
or defects, which are more favourable energetically than 
those at B sites. While the it band is strongly affected by 
such impurities, the a band is more robust, and there is 
little mixing of the two bands because of negligible inter- 
band scattering 0, ^| . This is consistent with accumu- 
lating experimental evidence suggesting that the a and ir 
bands are essentially in the ballistic and diffusive limit, 
respectively 0| (see also references in Refs. 0, Hi}). 
Nevertheless, so far in theoretical studies both of the two 
bands have been assumed to be either in the clean 0] 
or in the dirty limit |l7|. 

In this Letter, we examine theoretically the effects of 
induced superconductivity and impurities on the vortex- 
core structure in a two-band superconductor. Our work 
is motivated by the experimental investigation of the vor- 



tex state in MgB 2 using scanning tunnelling spectroscopy 
|lfj|. We present a novel formulation of coupled qua- 
siclassical Eilcnbcrger and Usadel equations to describe 
a multiband superconductor with one ballistic and one 
diffusive band with negligible interband scattering: the 
two bands are assumed to be coupled only by the pairing 
interaction. 

We apply our model to calculate numerically the lo- 
cal density of states (LDOS) and supercurrent density 
around an isolated vortex. We examine in detail the in- 
triguing spatial variation of these quantities and the order 
parameter in the two bands. A particularly interesting 
result emerging from our studies is the possibility of ad- 
ditional bound states near the gap edge in the a band. 

Our model is based on the equilibrium quasiclassical 
theory of superconductivity, where the physical informa- 
tion is contained in the Green function, or propagator, 
g(e, pfq, R). Here e is the quasiparticle energy measured 
from the chemical potential, pFa the quasiparticle mo- 
mentum on the Fermi surface of band a S {a, 7r}, and 
R is the spatial coordinate (the hat refers to the 2x2 
matrix structure of the propagator in the particle-hole 
space). In the clean a band, g a (e, Pf ct > R) satisfies the 
Eilcnberger equation 0] 



ef 3 - A CT , g a + vw ■ Vg a = 0, 



(1) 



where vp a is the Fermi velocity and A CT the (spatially 
varying) order parameter in the a band. Throughout 
this work, we ignore the external magnetic field (this is 
justified because MgB 2 is a strong type-II superconduc- 
tor). The coherence length in the a band is defined as 
6t = v f<t /2ttT c . 

In the presence of strong impurity scattering, the 
Green function has no momentum dependence and the 
Eilenberger equation reduces to the Usadel equation 0] . 
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We assume this to be the appropriate description for the 
dirty tt band, and take the propagator g^(e, R) to satisfy 



ef 3 — Aw , g-x 



-V • = 6. (2) 

TT 



The diffusion constant D defines the 7r-band coherence 
length as = ^D/2nT c . Additionally, both propaga- 
tors are normalized according to cfe = g% = —tt 2 1. 

We assume that the quasiparticles in different bands 
are coupled only through the pairing interaction, neglect- 
ing interband scattering by the impurities. The gap equa- 
tions for the multiband system are given as 



A a (R) = J2 VapN F pFp{R), 



(3) 



where a, (5 E {a, tt}, A a — fi Re A Q — T2 Im A a , the cou- 
pling matrix V a p determines the strength of the pairing 
interaction, N F p is the Fermi-surface density of states on 
band /3, and 



Yl (/ CT (e,PF<T,R)) PF „tanh^ 2r; 



J~ x(R) = 

J —e c 

^r(R) = / ' ^/-( £ ' R ) tailh (2^) 



(4) 



Here f a is the upper off-diagonal (f,2) element of the 
matrix propagator g a , {■■■) PFa denotes averaging over 
the a band Fermi surface, and e c is a cut-off energy. 

We solve the system of equations (0 (@J numerically. 
The normalization condition is taken into account with 
the Riccati parameterisation for the Green functions [2JJ, 
Elf . After self-consistency has been achieved for the order 
parameter, the (for the a band angle-resolved) LDOS in 
each band can be calculated as 



N a {e,p Fa ,R)/N Fa 
N w (e,R)/N F „ 



Im g a (e,p Fa ,R)/TT, 

Im gTr(e,TL)/w, (5) 



where g a is the upper diagonal (1,1) element of g a . 

The current density around the vortex has contribu- 
tions from both the tt band and the a band. The corre- 
sponding expressions are (e = — |e| is the electron charge) 



■U(R) 

2eN Fa 

.UR) 

2eN F ^ 



— (vpn-Im q,r) n ^ tanh ( — ^ , 
, 2tt x w/pf* \2TJ 



Throughout this work, we focus on the effects of purely 
induced superconductivity in the tt band. Diagonaliza- 
tion of the coupling matrix in Eq. decouples the gap 
equations. The larger of the two eigenvalues of the matrix 
VapNpp, denoted by A^ ^, determines the critical tem- 
perature T c and can be eliminated together with e c . The 
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FIG. 1: Magnitude of the order parameter, |A CTj ^(a;)|, in 
the a- and tt band, as a function of coordinate i on a path 
through the vortex center, at T = 0.1T C : (a) for different 
ratios £tt/£ct and fixed strength of the induced 7r-band gap, 
parameterized by the mixing ratio (see text), p = 0.3; (b) for 
different p and fixed ^l^o — 5. A weak effective repulsion 
due to Coulomb interaction was assumed in the subdominant 
A' 1 ' pairing channel, parameterized by A = —0.1 (see text). 



smaller eigenvalue A*- 1 ) is parameterized by the cut-off in- 
dependent combination A = A (0) AW /{X^ - A^). The 
pairing interactions have been calculated with ab initio 
methods, employing an electron-phonon coupling model 
0.0 ; from these studies, taking Coulomb repulsion into 
account |9(, we estimate A < 0.3 for MgB2. We present 
results for A = —0.1 (implying a weak repulsion in the 
subdominant A' 1 ' channel). The ratio of the magnitude 
of the bulk gaps p = | A^ ulk |/| Ajj ulk | near T c parameter- 
izes the strength of the induced superconductivity in the 
tt band; experimentally p ~ 0.3 in MgB2- For simplicity, 
we have set N Fa = N FlT in our calculations; the densi- 
ties of states of the two bands are indeed of comparable 
size H, 0] ■ A cylindrical Fermi surface was used for the 
two-dimensional a band. For the ratio of the coherence 
lengths in the two bands we present results for £„/£ a =l, 
3, and 5. An estimate from experiments [I0I.I22I gives for 
MgB2 a value between 1 and 3. 

In Fig. ^ we present the order-parameter magnitudes 
for each band as a function of coordinate x along a path 
through the vortex center. In Fig.n(a) we show the order 
parameter variation at T = 0.1T C for several coherence- 
length ratios £n-/£er an d § a P ratio p = 0.3. Surprisingly, 
we find that an increase of £„• results in an increase of 
the recovery length of the order parameter (the charac- 
teristic length over which the order parameter changes 
from zero to the bulk value) in the a band, while the 
recovery length in the tt band is barely affected. Thus, 
the recovery lengths in the two bands can differ consid- 
erably, even though the superconductivity in the tt band 
is induced by the a band. In Fig. ^ (b), the gap ratio 
p is varied for fixed ^/Ccr = 5. One can see that, apart 
from the well-known increase of the bulk A a /T c ratio 
with increasing bulk A v , the recovery length in the a 
band increases considerably with increasing p. The tt- 
band recovery length is also enhanced together with that 
of the a band if p > 0.3, as clearly seen for p — 0.5. For 
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FIG. 2: (a) Local density of states (LDOS) for the 7r band, 
as a function of energy e for different distances r from the 
vortex center. The spectra are shifted in vertical direction for 
convenience, (b) LDOS at the chemical potential, e — 0, as a 
function of r for different p and fixed T = 0.5T C , £,w/£,a = 3. 



p = 0.1, A OT near the vortex core is enhanced from its 
bulk value. 

To understand these effects we have performed calcu- 
lations for different values of A that characterizes the 
Coulomb repulsion in the ir band. We have found that 
for A > the ratio |A W |/|A CT | is reduced in the core re- 
gion with respect to its bulk value, however is enhanced 
for A < 0. This leads to a renormalization of the recovery 
lengths of the order parameters in both bands. 

In Fig.Elwe show the spectral properties of the it band. 
The LDOS, shown in Fig. 0(a), is flat at the vortex center 
(r = 0), in agreement with the experiment of Ref. |l0| . 
Outside the vortex core the BCS density of states is re- 
covered. The decay of the zero-bias LDOS as a function 
of radial coordinate is shown in Fig. |21 (b) , and compared 
with that of the a band [obtained from the data in Fig. 
(a)]. The decay length of the zero-bias DOS is clearly 
different for the two bands. For the it band it is given by 
t;^ y/ A a / An , and thus dominated by the parameter p. In 
the a band the length scale of the decay is £„, and thus 
shorter than that in the ir band. The existence of two 
apparent length scales in the LDOS was also reported in 
the case of two clean bands 0] and two dirty bands ^} • 

In Fig. [3] we show the vortex-core spectra in the a 
band. The LDOS, A^(e,R) = (iV CT (e, p F<T , R)) PF<T , as a 
function of energy for different distances from the vor- 
tex center is plotted in Fig. (a), showing the well- 
known Caroli-de Gennes-Matricon bound-state bands at 
low energies. The new feature of our model is the ad- 
ditional bound states in the vortex-core region near the 
gap edges, that are clearly visible in Fig. (a). This is 
in strong contrast with the case of a single-band super- 
conductor, where the spectrum near the vortex center is 
suppressed at the gap edges, showing neither a coherence 
peak nor additional bound states. We note that the self- 
consistency of the order-parameter profile is essential for 
discussing the presence of the bound states at the gap 
edges. We illustrate the development of these additional 
bound states in terms of the spectrum at the vortex cen- 
ter. We find that for a given £tt/£ct, the bound state exists 
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FIG. 3: (a) LDOS for the a band as a function of energy e 
for different distances r from the vortex center, (b) The de- 
velopment of an extra bound state near the gap edge in the a 
band for different parameter combinations. The bound state 
develops for sufficiently large mixing ratio p and coherence 
length ratio ^/Cu- (c) Angle-resolved LDOS for quasiparti- 
cles with momentum along the x axis at positions along the y 
axis in the a band, showing the bound states as a function of 
y. The dispersion of the bound states as a function of y is also 
shown in (d) for p — 0.5, ^tt/Cct = 3 and T — 0.5T C . Spectra 
in (a)-(c) are shifted in vertical direction for convenience. 



for p larger than a certain critical value, as seen in the 
left panel of Fig. 0(b). For a given p, on the other hand, 
the bound state develops if £ w /£cr exceeds a critical value, 
which is between 1 and 3 for p = 0.3, as can be inferred 
from the right panel of Fig. (b) . The bound states are, 
e.g., clearly resolved for p = 0.3 and £tt/£<t = 3, values 
appropriate for MgB2. We have found a similar bound 
state spectrum also for A > 0. The bound states move 
with the gap edge as a function of temperature. 

The bound-state spectrum is most clearly discussed 
in terms of the angle-resolved spectra, shown in Fig. 
(c) . Here the spectrum of quasiparticles moving in the x 
direction is shown as a function of position along the y 
axis. The position of the bound states as a function of 
y is shown in Fig. (d). The main bound state branch 
crosses the chemical potential in the vortex center. The 
additional bound states are seen near the gap edge, and 
show a weak dispersion. In fact, a close inspection reveals 
that there are two additional branches; however, only one 
of them is present at the vortex center. 

Finally, in Fig. 0] we show contributions from the two 
bands to the supercurrent density around the vortex 
separately. We show the results for p — 0.3, and (a) 

W&r = i) ( b ) £tt/£o- = 3 ) and ( c ) &t/6t = 5 > at 

T/T c = 0.1, 0.3 and 0.5. The current density from the 
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FIG. 4: Contribution to the current density as a function 
of distance from the vortex center, shown separately for the 
7T band and the a band for fixed mixing ratio p — 0.3 and 
different temperatures T; (a) for £,t/£ct = 1, (b) £tt/£ct = 3, 
and (c) ^tt/Co- = 5. In (d) the cr-band contribution to the total 
current density is shown as a function of radial coordinate, for 
T = 0.1T C , p = 0.3, and different &r/f CT . 



a band is enhanced for low temperatures near the vor- 
tex center, and the maximum approaches the center for 
T — > 0. This is due to the well-known Kramer-Pesch 
effect for the clean a band. The current density due 
to the induced superconductivity in the ir band is also 
enhanced by decreasing temperature, but the maximum 
does not approach the vortex center as in the a band. 
With increasing r, the contribution of the a band is re- 
duced, and the contribution of the 7r band becomes con- 
siderable when £,r/£ CT > 1. At the same time, \ a shows 
temperature dependence only in the core area, whereas 
Jtt is temperature dependent also far outside the core. 
To discuss this effect, we note that at a large distance r 
from the vortex center (but small compared to the Lon- 
don penetration depth), the current-density magnitudes 
are approximately given by j a (r) ~ eNp (T v Fa ./2r and 
jn(r) ~ eNp^TrDlA^l/r. The temperature dependence 
of in is thus dominated by that of |A T |. 

Another important observation is that, for the param- 
eters appropriate for MgB2 , the contribution of j n is con- 
siderable outside the vortex core. In fact, already for a 
moderate ratio ^/Co- = 3 the cr-band contribution is re- 
stricted to the region very close to the vortex center and 
is negligible outside the core for p > 0.3. To understand 
this effect, shown in Fig. 0(d), we note that the current- 
density ratio far away from the vortex center approaches 



lim 



jjf) = 2Dtt\A^\ _ \A„\ 



dominates when (£ n /£<r) 2 > (|A W |/T C ) 1 , which is a case 
relevant for MgB2. 

In conclusion, we have formulated a model for coupled 
ballistic and diffusive bands in terms of coupled Eilcn- 
berger and Usadel equations. We have studied the ef- 
fects of induced superconductivity in the 'weak' diffusive 
band on the order parameter, the current density and 
the spectral properties of the 'strong' ballistic band. We 
have found that (a) the recovery lengths of the order pa- 
rameters in the two bands are renormalized by Coulomb 
interactions; (b) the vortex core spectrum in the a band 
shows additional bound states at the gap edges; and (c) 
the current density is dominated in the vortex core by 
the cr-band contribution, and outside the vortex core the 
7r band contribution is substantial, or even dominating, 
for parameters appropriate for MgB2. Our predictions 
concerning the vortex-core spectrum of the a band can 
be tested in future tunnelling experiments. 
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